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The  Coupled  Damping  Coefficients 
of  a  Symmetric  Ship 


By  R.  Tim  man1  and  J.  N.  Newman1 


A  study  h  mod*  of  a  floating  or  tubmergsd  body  with  longitudinal  and  transverse  sym- 
metry^whkh  is  moving  with  constant  forward  speed  and  performing  small  oscillations. 
The  analysis  is  quite  general  in  the  sense  that  the  shape  of  the  body  and  the  nature  of  the 
oscillations  are  unspecified,  but  it  is  assumed  that  the  linearized  free-turface  condition 
holds.  With  this  assumption  the  oscillatory  velocity  potential  is  found  in  terms  of  on 
unknown  Green's  function,  the  existence  of  which  is  also  assumed.  This  potential  is  then 


used  to  show  the  symmetry  properties  of  the 


A  controversy  has  arisen  in  ship-motion  theory  re¬ 
garding  the  cross-coupling  damping  coefficients  of  a 
pitching  and  heaving  ship.  If  the  ship  is  symmetrical 
fore  and  aft  and  if  it  has  no  forward  speed,  then  from 
symmetry  the  cross-coupling  moment  due  to  heave  and 
the  force  due  to  pitch  must  both  be  zero,  at  least  in  the 
linearized  solution  of  the  problem.  However,  if  the  ship 
is  moving  with  forward  speed,  an  asymmetry  is  intro¬ 
duced  and  cross-coupling  results.  Haskind  [1]’  has  em¬ 
ployed  thin-ship  theory  to  show  that  for  a  symmetric 
ship  with  constant  forward  speed,  the  two  cross-coupling 
damping  coefficients  for  pitch  and  heave  are  equal  in 
magnitude  and  opposite  in  sign.  However,  this  con¬ 
clusion  has  received  criticism  in  several  papers  and  a  dis- 
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cross-coupling  damping  coefficients. 

pute  has  arisen,  which  is  reviewed  in  the  survey  of  Vossers 

[2]- 

An  analysis  of  particular  relevance  to  this  discussion  is 
that  of  Havelock  [3],  which  considers  the  case  of  a  float¬ 
ing  spheroid  with  a  rigid  free-surface  condition.  Have¬ 
lock  assumes  that  the  spheroid  is  pitching  and  heaving 
and  that  there  is  a  constant  forward  speed  or,  equiva¬ 
lently,  that  there  is  a  uniform  flow  of  the  stream.  With 
this  model  it  is  found  that  the  two  cross-coupling  co¬ 
efficients  are  of  unequal  magnitude,  and  Havelock  con¬ 
cludes  that  equality  of  the  cross-coupling  is  a  consequence 
of  the  thin-ship  approximation.  However,  the  sum  of 
the  two  coefficients  can  be  expressed  [4]  in  terms  of 
energy  radiated  in  outgoing  surface  waves  and,  if  a  rigid 
free  surface  is  assumed,  there  can  be  no  waves  and  there¬ 
fore  no  energy  radiation.  This  reasoning  leads  to  a  con¬ 
tradiction  with  Havelock’s  result,  for  if  the  sum  of  the 
two  cross-coupling  coefficients  is  zero,  then  they  must  be 
equal  and  opposite. 

The  source  of  this  discrepancy  lies  in  the  fact  that  in 
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damping  coefficient* 
forward  velocity 
equation  of  the  body  surface 
co*(n,  z)  -  horizontal  direc¬ 
tion  coeine 

COW n,  y)  -  transverse  direc¬ 
tion  cosine 

coe(a,  s)  -  vertical  direction 
coeine 

V  co*(n, «)  -*eoe(n,y) 
a  oo*(n,  z)  -  *  co*  (n,  a) 
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t  co My)  -  y  co*(n,  z) 
Green's  function 
gravitational  acceleration 
unit  vectors 

unit  normal  into  body  surf&oe' 

fluid  preeture 

time 

velocity  vector  of  steady  flow 
field 

Cartesian  co-ordinates  fixed 
in  space 


(*',  y',t')  -  Cartesian  co-ordinate*  fixed 
in  body 

a  -  displacement  vector  of  a  point 
on  body 

fs  -  oscillation  amplitudes  of  body 
f,  S,  f  -  dummy  oo-ordinate*  corre¬ 
sponding  to  (z,  y,  t) 

M)  m  velocity  potential  of  oscilla¬ 
tory  flow 

w  -  circular  frequency  of  oscilla¬ 
tions 


Havelock’s  analysis  the  boundary  condition  is  satisfied 
by  taking  the  oscillatory  normal  velocity  on  the  spheroid 
and  equating  this  to  the  normal  velocity  of  the  fluid  on 
the  mean  potition  of  the  spheroid  in  space.  In  fact,  the 
oscillatory  disturbance  is  a  small  perturbation  of 
the  steady  flow  field,  and  the  boundary  condition  on  the 
spheroid  must  be  satisfied  on  the  exact  oscillating  surface 
of  the  body,  or  else  expanded  to  the  mean  surface  in  a  sys¬ 
tematic  manner  so  as  to  include  the  oscillatory  flow  in¬ 
duced  on  the  body  surface  by  its  change  of  position  in  the 
steady-state  field.  That  is,  the  oscillations  of  a  ship  in  a 
moving  stream  give  rise  to  a  small  disturbance  of  the 
steady  flow  field,  and  various  second-order  effects  enter 
into  the  unsteady  problem  as  a  result  of  the  lower-order 
steady  field.  It  can  be  shown  by  an  extension  of  Have¬ 
lock’s  analysis  that  if  the  boundary  condition  is  satisfied 
on  the  exact  surface  of  the  spheroid,  then  the  cross-cou¬ 
pling  damping  coefficients  between  pitch  and  heave  are  in 
fact  equal  in  magnitude  and  opposite  in  sign.  It  thus 
seems  plausible  that  this  equivalence  holds  for  any  sym¬ 
metrical  ship  or  body,  irrespective  of  the  thin-ship  as¬ 
sumption. 

In  order  to  study  this  question  more  generally,  the 
present  paper  treats  the  problem  of  an  arbitrary  floating 
or  submerged  body  with  longitudinal  and  transverse 
symmetry,  which  is  moving  with  constant  forward  speed 
and  oscillating  sinusoidally  in  any  of  the  six  degrees  of 
freedom.  The  only  significant  assumptions  arc  that  the 
problem  is  linear,  in  the  sense  that  the  oscillations  are 
small  and  that  the  disturbance  of  the  free  surface  due  to 
the  forward  motion  is  also  small.  To  be  physically 
realistic,  the  latter  assumption  implies  that  the  body  is 
thin,  slender,  or  deeply  submerged,  or  a  combination  of 
these,  but  the  analysis  and  conclusions  are  equally  valid 
for  the  case  of  a  nonslender  body  with  a  rigid  free-surface 
condition.4  (“Shallow"  ships,  with  small  draft  and 
finite  beam,  are  not  included  in  the  present  work.)  With 
these  two  basic  assumptions  we  show  that  the  sum  of  the 
two  complementary  cross-coupling  damping  coefficients 
is  zero  for  all  pairs  of  modes  of  oscillation  except  for  the 
coupling  of  surge  with  pitch  and  roll  with  sway.  Further¬ 
more,  fifteen  of  the  thirty  cross-coupling  coefficients  are 
shown  to  be  zero.  The  same  conclusions  have  been  ob¬ 
tained  for  a  thin  ship  by  Hanaoka  [6], 

As  usual  we  assume  irrotational  incompressible  flow 
and  formulate  the  problem  in  terms  of  the  velocity  po¬ 
tential.  The  potential  problem  is  solved  in  terms  of  a 
Green’s  function  which  is  not  explicitly  known,  but  the 
existence  of  this  function  seems  physically  plausible  and 
can  probably  be  proved  by  recourse  to  the  theory  of 
Fredholm  integral  equations.  The  reciprocity  properties 
of  this  Green’s  function  are  then  established  and  the 


symmetry  properties  of  the  damping  coefficients  follow 
directly.  It  should  be  emphasized  that  we  do  not  ex¬ 
plicitly  solve  either  the  steady  or  unsteady  potential 
problems. 

The  proper  representation  of  the  oscillating  potential, 
equation  (9),  is  particularly  interesting,  since  it  demon¬ 
strates  the  effect  of  satisfying  the  boundary  condition  on 
the  exact  (oscillating)  surface  of  the  body.  The  final 
results  for  the  dumping  coefficients  are  shown  in  a 
matrix,  Table  1. 

Both  the  analysis  and  final  results  are  analogous  to 
reciprocity  studies  in  aerodynamics  [7,  8],  In  fact  this 
analogy  was  the  original  motivation  for  suspecting  that 
the  equivalence  of  the  pitch  and  heave  cross-coupling  co¬ 
efficients  did  not  depend  on  the  thin-ship  assumption. 

The  Boundary-Value  Problem 

Let  (x,  y,  z)  be  a  Cartesian  co-ordinate  system,  mov¬ 
ing  through  the  fluid  with  constant  velocity  c,  with  z  ver¬ 
tically  upward  and  x  in  the  direction  of  forward  mo¬ 
tion.  In  addition  we  shall  employ  an  oscillatory  co-ordi¬ 
nate  system  *'  =  *  —  of’"*  where  a  is  an  infinitesimal 
vector,  which  may  depend  on  and  the  real  part  is  to 
be  taken  in  expressions  involving  e The  ji'  co-ordi¬ 
nates  are  fixed  with  respect  to  a  body  which  is  defined  by 
the  equation  F(x',  y',  z')  -  0.  The  velocity  of  the  fluid 
is  represented  by  the  vector 

»(»)  +  e'"'v4>(*) 

Thus  v  is  the  steady  velocity  field  due  to  the  forward 
motion  of  the  body,  in  the  presence  of  the  free  surface, 
and  *(x)  is  the  potential  of  the  oscillating  velocity  vector. 
The  function  4>  must  satisfy  Laplace’s  equation,  the 
linearized  free-surface  condition  [4] 

-9  ~  +  2m*  ^  -  c*  ^  =  0  (1) 

Oz  Ox  ox* 


on  the  undisturbed  free  surface  z  =  0,  and  a  suitable 
radiation  condition  at  infinity. 

The  boundary  condition  on  the  body  is 
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‘The  distinction  between  thin  and  elender  bodies  is  important. 
In  both  cases  the  beam  is  small  compared  to  the  length  but  the 
thin  ship  has  a  small  beam-draft  ratio  as  well,  whereas  the  slender 
ship  has  beam  and  draft  of  the  same  order  of  magnitude.  By 
deeply  rubmerged  we  imply  that  the  depth  of  submergence  is  suffi¬ 
ciently  large  that  the  waves  will  be  small,  but  not  so  large  that 
the  free  surface  may  be  neglected  and  the  fluid  considered  a a  in¬ 
finite. 
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This  condition  holds  on  the  actual  surface  of  the  body 
where,  from  Taylor’s  theorem 

»  =  (v).„.  +  [(e-v)v]„,.,«'“'  +  C(a’) 

The  subscript  “mean”  denotes  that  the  fu  Miction  is  to  be 
evaluated  on  the  mean  position  of  the  fcody,  or  with 
x  —  x'.  Thus  we  obtain  the  boundary  eoncMition 

0  =  -W(«V,.F)  +  [v(x')  + 

0{at) 

The  steady-state  term  gives  the  boundary  condition  for 

v, 

v(x')-V,'F  =  0,  (2) 

and  the  terms  of  first  order  in  the  small  oscillatory  func¬ 
tions  give  the  boundary  condition  for  4, 

X4-V.F  =  iua-v.F  -  |(«  T)v]m.«  rr  -t 

+  M—  ■  [l  V.pj  +  i  (tr  Vrf'j 

+fc(s  •’-')] 

or 


^(x)  =  ^  //  [G(m-  °^/(0 

-  ♦(0^G(*,0]d««  (5) 

where  x  is  a  point  on  the  mean  body  surface  and  the  in¬ 
tegration  is  over  this  surface,  the  undisturbed  free  surface 
*  =  0,  and  a  closing  surface  at  infinity.  The  Green’s 
function  is  any  harmonic  function  of  x  and  (  which  is 
singular  like  the  potential  for  a  source  at  the  point  x  = 
We  shall  assume  the  existence  of  a  Green’s  function  which 
satisfies  the  free-surface  condition,  the  same  radiation 
condition  is  4,  and  the  condition 

—  =  0  on  the  body  F (?)  =  0 

As  stated  in  the  introduction,  the  existence  of  this 
Green’s  function  can  probably  be  proven  from  the  theory 
of  Fredholm  integral  equations,  and  furthermore  seems 
physically  plausible  since  this  function  can  be  visualized 
as  the  potential  due  to  a  realistic  fluid  flow;  i.e.,  the 
disturbance  caused  by  a  small  pulsating  sphere  in  the 
presence  of  the  body  and  the  free  surface. 

Before  proceeding  further  we  must  establish  the  re¬ 
ciprocal  properties  of  this  Green’s  function.  Assume  that 
G *  and  G~  are  two  Green’s  functions  satisfying  the  con¬ 
dition  dG/6n,  =  0  on  the  body  and  satisfying  the  free- 
surface  conditions 


V4V.-F  *=  iua-V.-F 


+  [(*■»)«  -  (*’  V)v]  WTWIU'  XsrF  (3) 

In  equation  (3)  and  hereafter,  all  of  the  t-«rms  are  small, 
of  the  same  order  as  a  or  4-  Thus  to  tbiis  order  of  ap¬ 
proximation  it  is  no  longer  necessary  to  distinguish  be¬ 
tween  the  actual  position  of  the  body  and  its  mean 
position,  or  between  the  co-ordinates  i  and  x'. 

We  now  use  the  vector  identity* 

(*'?)«  ~  («•  V)v  =  V  X  («  X  v)  -  xT  i  -f  vV-« 

and  since  v  y  =  0  (from  incompressibility)  and  y-vF 
=  0  [from  equation  (2)],  we  find  that 

V4VF  =  iua-vF  +  [v  X  (o  X  v)]-T^  on  F  =  0 


+  o,«G±  =F  2iwc^ 
df  o( 


c* 


d«G± 

df* 


«■  0  on  f  =  0.  (6) 


These  two  functions  correspond  physically  to  the  velocity 
potentials  of  an  oscillating  source  at  the  point  x  =  {,  in 
the  presence  of  the  body  and  the  free  surface.  The  func¬ 
tion  G*  corresponds  to  the  case  where  there  is  a  free- 
stream  velocity  c  in  the  —  x-direction  and  the  function 
G~  to  a  flow  with  velocity  c  in  the  -fx-direction.  Be¬ 
cause  of  this  difference  the  sign  of  the  third  term  in  (6) 
must  differ  for  the  two  cases. 

From  Green’s  theorem,  with  x  and  y  two  different 
vectors, 


Since  vF  is  a  vector  normal  to  the  Ixxly  surface,  it 
follows  that  the  boundary  condition  for  4  on  the  body 
may  be  written  as 

r  =  [«««  4-  V  X  (a  X  *)]  •«  (4) 

On 

Thus  the  effect  of  the  steady  flow  is  "to  increase  the 
normal  oscillating  velocity  by  v  X  (e  X  y)  • 


G+(x,  y)  -  G-(y,  x) 

-  C-(y,  0vG+(x,  0]dF( 

■  s//  0  s,  °  -  °'w’  ° 


Green's  Theorem  and  tha  Green's  Function 

Now  we  employ  Green’s  theorem 


*Cf,  [9],  equation  (1.4.13). 


where  the  surface  integral  is  over  the  body,  the  undis¬ 
turbed  free  surface,  and  a  closure  at  infinity.  The  inte- 
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gral  on  the  body  vanishes  since  dG^/dn <  -  0  and  the 
closure  at  infinity  vanishes  from  the  radiation  condition.' 
The  integral  over  the  free  surface  is  equal  to 

+  G~  (uxO+  +  c'~  d(d, , 

-  2-^/0*  H'-'-f) 

+  G~  (uocG*  +  c‘^)]rfl 

where  the  line  integral  is  over  the  boundary  or  boundaries 
or  the  free  surface,  or  the  intersection  of  the  free  surface 
with  the  closure  at  infinity  and  (if  any)  with  the)  body. 
From  the  radiation  condition  the  integral  over  the, boun¬ 
dary  at  infinity  vanishes,  and  if  the  body  is  submerged 
there  is  no  further  boundary,  with  the  result  that 

G+(z,  V)  -  o-(v,  s)  -  0 
or 

G+(»,  «  -  <?-«,  *) 

For  a  floating  body,  the  beam  must  be  small  (i.e.,  the 
body  is  either  thin  or  slender)  and  thus  the  line  integral 
around  the  waterline  is  of  order 


on  the  undisturbed  free  surface,  and  thus  that  satisfies 
the  free-surface  condition  (1).  Furthermore,  the  surface 
integral  in  (5)  may  be  treated  in  exactly  the  same  manner 
as  we  did  in  establishing  the  reciprocity  relation  between 
the  Green’s  functions.  It  follows  that 

*(*)  -  ^yy  g+(*.  o  ^  *«)<«?*  (8) 

where  the  integration  is  only  over  the  body. 

Substituting  (4)  in  (8)  we  obtain 

*  ”  ~r  jy  G(x,  ()  [tua  +  V  X  («  X  v)  [ndSi 

This  is  equal  to 

d  m  jy  {twaG  +  V  X  ((«  X  v)G) 

-  (VG)  X  («  X  v) }  ■  ndS 

but  from  Stokes’  theorem 

SSv  X  [(«  X  v)G]rfS  =  »  X  v)G  <f| 

where  the  line  integral  is  again  over  the  intersection,  if 
any,  of  the  body  with  the  undisturbed  free  surface. 
Once  again  we  invoke  the  linearized  free-surface  con¬ 
dition  ;  if  the  waves  are  small,  then  v  on  the  free  surface  is, 
to  first  order,  tangent  to  the  undisturbed  plane  of  the 
free  surface.  Since  v  is  also  tangent  to  the  body  surface, 
it  is  tangent  to  the  intersection  l,and  thus  («  X  v)  dl  -  0, 
and  the  line  integral  vanishes.  Thus  we  find  that 

♦  -  j^yy  {tu«G  -  VG  X  («  X  v))ndS 

or,  since 

[VG  X  (•  X  »)]•«  =  [(VGv)«  -  («  VG)v]  ii 

-  (wG)(«  n) 

it  follows  that 


/N 

which  is  of  the  same  order  as  the  beam.  Thus  it  is  con¬ 
sistent  with  the  linearized  free-surface  condition  that,  in 
all  cases, 

G+(k,  0  =  G~«,  a)  (7) 

That  is,  the  Green’s  function  is  reciprocal  if  the  direction 
of  the  streaming  flow  is  reversed.  This  property  is  well 
known  for  the  Green’s  function  which  does  not  satisfy  a 
boundary  condition  on  the  body. 

We  now  return  to  the  construction  of  the  velocity 
potential  from  equation  (5)  substituting  the  Green’s 
function  G+(x,  0-  From  the  reciprocal  property  (7)  it 
follows  that 


dC+  ,  dG+  ,d«G+ 

—a  — - h  (• )*G+  +  2twc  — - c«  - 

dt  dx 


dz» 


•  This  consequence  of  the  radiation  condition  is  not  physically 
obvious.  Some  discussion  of  this  point  will  be  found  in  reference 

(10],  page  458. 


♦  =  h  JJ (*»°  -  wiG)(„.n)dS,  (9) 

where  the  integral  is  over  the  body  surface.  Thus  the 
effect  oj  the  steady  velocity  field  on  the  unsteady  potential  is 
expressed  by  the  factor  v-vG.  Physically  this  can  be 
thought  of  as  a  dipole  distribution  in  the  direction  of  v 
(and  thus  tangent  to  the  body  surface)  and  of  strength 
equal  to  the  normal  displacement  («-n)  times  the  mag¬ 
nitude  of  v.  The  same  result  has  been  derived  for  a  thin 
ship  [4]. 

The  Forces  ond  Moments 

Equation  (9)  holds  for  any  oscillatory  displacement 
vector  o,  and  is  therefore  not  restricted  to  rigid  body 
motions.  We  now  assume  that  the  body  is  rigid,  with 
six  degrees  of  freedom.  It  is  convenient  to  introduce  an 
indicial  notation,  where  we  denote  the  six  oscillatory 
velocities  by 

*'« M*  O'  -  1,  2,  ...  6) 
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These  are,  respectively,  surge,  sway,  heave,  roll,  pitch, 
and  yaw.  We  also  define  the  six  matrix  elements 

/)(«)  =  cos(n,  x) 

/.(*)  =  cos(n,  y) 


the  direction  of  the  forward  velocity  reversed.  To  dis¬ 
tinguish  between  these  two  eases  we  shall  denote  the 
matrix  appropriate  to  forward  motion  in  the  +  x-direc- 
tion  as  B„*  and  the  matrix  of  the  reverse  flow  by  B„~. 
Then 


/,(*)  =  cos(n,  z) 

/«(*)  =  y  cos  (n,  t)  —  t  cos  (n,  y) 

/%(*)  “  t  cos  (n,  x )  —  x  cos  («,  z) 

/»(«)  -  x  cos  (n,  y)  —  y  cos  (n,  x) 

Then  the  normal  displacement  a-n  at  the  point  x  is  given 
by 

«  n  =  £  fiA(*)  00) 

j'-i 

and  if  p  is  the  hydrodynamic  pressure,  the  hydrodynamic 
forces  and  moments  are  given  by  the  six  expressions 

F,  -  -  ffpf,(»)dS ,  («  =  1,  2,  ...  6)  (11) 

The  pressure  p  is,  from  the  linearized  form  of  Bernoulli’s 
equation, 

p  -  —  p[iu^e'“‘  +  fvV*)e'“‘  +  M  (*'*)|  (12) 

where  terms  of  second  order  in  the  oscillatory  potential  * 
are  neglected. 

The  integral  in  (11)  must  1*  evaluated  on  the  oscillat¬ 
ing  surface  of  the  ship,  and  thus  the  zero-order  term 
x(v  v)  must  be  expanded  to 

W(»'»)w>  = 

+  «-(rT)(M*-*)i»«  +  0(«’) 

Also  if  the  lx>dy  intersects  the  free  surface,  the  oscillatory 
change  in  the  surface  of  integration  must  be  included. 
However,  both  of  these  effects  are  in  phase  with  the  dis¬ 
placement  oe'"1  and  will  not  influence  the  damping  co¬ 
efficients. 

We  restrict  ourselves  then  to  the  damping  forces  and 
moments,  which  may  be  represented  by  the  imaginary 
part  of  the  integral  (11),  taken  over  the  mean  surface  of 
the  body.  let  the  damping  coefficients  be  represented 
by  the  matrix  B,„  where  the  first  index  denotes  the 
direction  of  the  force  and  the  second  index  the  velocity 
component  involved.  Thus,  for  example,  is  the  heave 
damping  force  due  to  pitching  oscillations.  The  thirty- 
six  coefficients  are  identified  by  the  matrix  in  Table  1. 
Combining  equations  (9-12)  and  taking  the  imaginary 
part,  we  obtain  the  expressions 

B„  =  ^  Im  Jj  /<(x)(to>  +  v(x)  •▼,) 


/,«)(*«  -  v«)-Tt )G+(x,  t)dS,dS,  (13) 


In  this  form  the  possibility  of  symmetry,  and  the  im¬ 
portance  of  the  dipole  distribution  — v  VjG  are  apparent. 

In  order  to  establish  the  symmetry  properties  of  the 
coefficients  we  consider  the  same  forces  and  moments  with 


B,,±  =  2t  lmJI 

ff  JA ©(*  -  vk(t)  r,)G±(«,  0 dS<dS.  (14) 

where  6,±(x,  0  arc  the  two  Green’s  functions  which  we 
introduces!  Ixdore,  and  which  possess  the  reciprocal 
property 

G+(x,  0  «  G-(«,  x) 

We  need  one  further  assumption  based  upon  the 
linearized  frec-surface  condition.  This  is  that,  to  first 
order, 

v+(>:).=  -v(x) 

on  the  body.  This  seems  consistent  with  the  assumption 
of  a  small  disturbance  on  the  free  surface,  for  if  the  waves 
are  small,  the  effect  of  the  free  surface  on  the  flow  at  the 
body  is  small,  and  thus  to  first  order  the  steady  velocity 
on  the  body  is  an  odd  function  of  the  stream  velocity  c,  as 
is  the  case  for  a  body  in  a  wave-free  field. 

Substituting  v~  ■=  — v+  and 

C~(x,  0  -  G+«,  x), 

it  follows  that 

B,r  =  /,<«)(<■#  ~  v+(*)  ▼.) 

or,  after  interchanging  the  integrals  and  the  variables  of 
integration  and  comparing  the  resulting  expression  with 
(14), 

B,r  =  B,S  (15) 

It  should  be  noted  that  up  to  this  point  in  the  analysis, 
no  assumption  has  been  made  regarding  the  symmetry  of 
the  body.  Thus  (15)  holds  for  asymmetric  bodies  and, 
in  particular,  it  follows  that  the  six  principal  damping  co¬ 
efficients  B„  are  independent  of  the  direction  of  forward 
motion.  Of  greater  practical  importance  is  the  fact  that 
if  there  is  no  forward  speed,  Bt,+  and  B„~  must  be  the 
same,  and  therefore  B„  =  B,,  when  c  =  0.  This  result 
has  also  been  obtained  by  Haskind  [11]  and  confirms 
the  frequent  argument  based  upon  strip  theory. 

Now  we  consider  the  physical  relations  between  B„+ 
and  B„~,  assuming  that  the  ship  is  symmetrical.  For 
example,  the  pitch  moment  of  a  symmetrical  body  due  to 
heave  and  heave  force  due  to  pitch  are  odd  functions  of 
the  forward  velocity,  since  changing  the  direction  of  the 
flow  is  equivalent  to  looking  at  the  body  from  the  op¬ 
posite  side  and  changing  the  sign  of  rotation  about  the 
pitch  axis.  Examination  of  all  of  the  coefficients  in  this 
manner  yields  the  following  conclusions: 
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Table  1  Matrix  of  Symmetry  Propartist  of  Damping  Crott-Coupling  Coefficients 


Surge 

Sway 

Heave 

Roll 

Pitch 

Yaw 

Surge  force 

i  -  i 

3  ~  2 

i  -  i 

3  "  4 

3  -  5 

i  -  a 

i  -  1 

Swav  force 

flu 

flu  -  0 

flu  -  -flu 

flu-0 

flu  —  fl.  1 

flu-0 

i  -  *2 

Heave  force 

«,,  -0 

flu 

flu-  0 

flu  —  fl« 

fl,. -0 

flu  —  —  flu 

«  —  3 

Roll  moment 

flu  ■  —  fl„ 

flu  -  0 

flu 

flu  -  0 

flu  —  —flu 

flu -0 

i  —  4 

Pitch  moment 

«„  -0 

flu  -  fl;. 

flu  -0 

fl.. 

fl.,  =  0 

flu  -  -flu 

i  -  5 

Yaw  moment 

«„  -  flu 

flu  -  0 

flu  —  —  flu 

flu  -  0 

«u 

flu  -  0 

i  -  « 

fl.,  -  0 

fl.;  -  -Hu 

flu  -  0 

flu  —  —  flu 

flu  -  0 

flu 

1  The  six  coefficients  B„  and  the  nine  coefficients 
flu,  flu,  flu,  fl«,  flu,  B«,  flu,  flu,  and  flu  arc  even  funet  ions 
of  the  forward  velocity  c. 

2  The  twelve  coefficients  flu,  flu,  flu,  flu,  B«,  fl»«, 
flu,  fl«,  flu,  flH,  fln,  and  flM  are  odd  functions  of  e. 

3  The  nine  coefficients  flu,  Bn,  fl»,  flu,  flu,  flu,  flu, 
flu,  and  flu  are  all  zero  by  symmetry;  i.c.,  there  are  no 
transverse  forces  or  moments  due  to  longitudinal  oscilla¬ 
tions. 

Combining  these  conclusions  with  equation  (15)  it 
follows  that: 

4  The  sum  B„  +  fl,,  =  0  for  all  cross-coupling  co¬ 
efficients  except  flu  +  flu  (su:gc  and  pitch)  and  B«  + 
fl«  (roll  and  sway),  where  B„  —  B,,  -  0. 

5  Half  of  the  cross-coupling  coefficients  are  zero 
(those  listed  under  statement  3  plus  their  complementary 
members).  Furthermore  this  conclusion  holds  for  an 
asymmetric  body  as  well,  since  statement  3  and  equation 
(15)  are  valid  without  the  assumption  of  longitudinal 
symmetry. 

Discussion  of  Results 

The  symmetry  properties  of  the  damping  coefficients 
for  a  symmetrical  ship  are  shown  in  Table  1.  For  the 
cross-coupling  coefficients  it  is  seen  that  B„  4-  fl,,  =  0 
except  for  coupling  between  pitch  and  surge,  where 
flu  =  flu,  and  for  coupling  Irctween  roll  and  sway,  where 
Bu  =  flu.  These  results  are  consistent  with  the  con¬ 
clusions  based  upon  thin-ship  theory  [1,  4,  5,  6],  but  the 
present  derivation  is  valid  for  all  thirty  cross-coupling 
damping  coefficients  and  docs  not  require  that  the  ship  be 
thin.  We  have,  however,  assumed  that  the  waves 
created  by  the  forward  motion  are  small,  and  this  implies 
that  the  body  is  either  thin,  slender,  or  deeply  submerged. 
Nevertheless  the  analysis  and  conclusions  also  hold  for 
Havelock’s  [3]  mathematics)  model  of  a  nonslender  body 
with  a  rigid  free  surface,  which  may  correspond  physically 
to  very  slow  forward  speed  and  a  low  frequency  of 
oscillations. 

There  does  not  appear  to  be  sufficient  experimental 
evidence  to  support  these  conclusions  completely,  but 
the  oscillator  experiments  of  Gerritsma  [12]  and  Golovato 
[13]  are  strongly  suggestive  of  the  equivalence  between 


the  pitch  and  heave  coupling  coefficients.  Gerritsma  has 
measured  t  he  coefficient  s  fl»  and  flM  for  a  Series  CO  model. 
Since  this  hull  is  not  longitudinally  symmetric,  the 
present  theory  is  not  strictly  valid.  Following  a  sugges¬ 
tion  of  Vossers  [2],  however,  we  may  separate  the  cross¬ 
coupling  coefficients  into  two  parts 

B„  -  fl,/  +  fl,/1 

where  B,/  is  the  value  of  B,,  at  zero  forward  speed  and 
fl,/1  is  the  difference  due  to  the  effects  of  forward  speed. 
From  equation  (15)  it  follows  that  (even  for  an  asym¬ 
metric  body) 

*«/  =  Btt 

If  the  ship  is  only  slightly  asymmetric,  the  effects  of  the 
asymmetry  on  the  coefficients  B,/1  will  be  small  and  thus, 
approximately, 

Bu"  -  -flu" 

The  experimental  results  in  Fig.  5  of  reference  [12]  do 
not  confirm  this  relation  exactly,  but  they  do  suggest, 
especially  at  the  higher  frequencies  and  speeds,  that  this 
is  a  meaningful  approximation  for  a  slightly  asymmetric 
ship. 

Experimental  measurements  of  Bu  and  flu  have  been 
made  for  a  symmetrical  hull  by  Golovato  (the  results  of 
the  pitch  experiments  are  unpublished),  who  found  that 
these  two  coefficients  were  approximately  equal  and 
opposite,  with  a  maximum  difference  of  about  10  per  cent 
over  a  fairly  wide  range  of  speeds  and  frequencies. 

With  regard  to  the  derivation  of  this  theory  an  im¬ 
portant  result  is  the  effect  of  the  steady  flow  field  on  the 
oscillatory  potential.  In  equation  (9)  this  effect  is  seen 
to  be  a  tangential  dipole  distribution,  equal  in  strength 
to  the  product  of  the  steady  velocity  and  the  normal  os¬ 
cillatory  displacement.  The  consideration  of  this  effect 
is  vital  to  the  present  analysis  as  it  is  directly  responsible 
for  the  symmetry  properties  of  the  cross-coupling  co¬ 
efficients.  Physically  this  implies  that  the  problem  of  an 
oscillating  body  in  a  moving  fluid  is  not  the  same,  even 
in  the  linearized  sense,  as  the  problem  of  a  fixed  body 
with  the  same  distribution  of  normal  velocity,  for  there  is 
an  additional  effect  from  moving  about  in  the  steady 
flow  field. 
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